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ABSTRACT 

The most common means of converting an observed CO line intensity into a molecular 
gas mass requires the use of a conversion factor (Xqq). While in the Milky Way this quantity 
does not appear to vary significantly, there is good reason to believe that Xqo will depend 
on the larger-scale galactic environment. With sensitive instruments pushing detections to 
increasingly high redshift, characterising Xqq as a function of physical conditions is crucial 
to our understanding of galaxy evolution. Utilising numerical models, we investigate how 
varying metallicities, gas temperatures and velocity dispersions in galaxies impact the way 
CO line emission traces the underlying H2 gas mass, and under what circumstances Xqq may 
differ from the Galactic mean value. We find that, due to the combined effects of increased 
gas temperature and velocity dispersion, Xco is depressed below the Galactic mean in high 
surface density environments such as ULIRGs. In contrast, in low metallicity environments, 
Xco tends to be higher than in the Milky Way, due to photodissociation of CO in metal- 
poor clouds. At higher redshifts, gas-rich discs may have gravitationally unstable clumps that 
are warm (due to increased star formation) and have elevated velocity dispersions. These discs 
tend to have Xco values ranging between present-epoch gas-rich mergers and quiescent discs 
at low-z. This model shows that on average, mergers do have lower Xco values than disc 
galaxies, though there is significant overlap. Xco varies smoothly with the local conditions 
within a galaxy, and is not a function of global galaxy morphology. We combine our results to 
provide a general fitting formula for Xco as a function of CO line intensity and metallicity. 
We show that replacing the traditional approach of using one constant Xco for starbursts and 
another for discs with our best-fit function produces star formation laws that are continuous 
rather than bimodal, and that have significantly reduced scatter. 

Key words: ISM:clouds-ISM:molecules-galaxies:interactions-galaxies:ISM- 
galaxies;starburst-galaxies:star formation 



1 INTRODUCTION 

As the building block of stars, H2 is arguably the most important 
molecule in astrophysics. Ironically, however, it is also one of the 
more observationally elusive. With no permanent dipole moment, 
H2 is best directly detected via its first quadrupole line. This line 
lies at ^ 500 K above ground, significantly above the ~ 10 K 
typical of the cold molecular interstellar medium (ISM), and in a 
spectral region with relatively low atmospheric transmission. As a 
result, giant molecular clouds (GMCs) are often studied via tracer 



molecules. The ground-state rotational transition of carbon monox- 
ide (^^co (1=1-0), hereafter CO) is one of the most common tracers 
of H2 in GMC s owing to its rel atively high abundance (~ 10~''/H2 
in the Galaxy: iLee et alj|l996h , the high atmospheric transmission 
at ~ 3 mm where the J=l-0 line lies, and the low temperatures and 
densities required for CO excitation (~ 5 K, ~ 10^ — 10'^ cm^'^). 
However, using CO to trace H2 does not come without uncertainty. 
At the basis of the interpretation of CO observations is the con- 
version between CO spectral line intensity and H2 gas mass, the 
so-called CO-H2 conversion factor. 

The CO-H2 conversion factor is defined as 
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where A'^h2 is the molecular gas column density, and Wco is the 
velocity-integrated CO line intensity (measured in K-km s"'^). Al- 
ternatively, the conversion factor can be defined as the ratio of the 
molecular gas mass and CO line luminosity: 
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Xco and aco are easily related via 

Xco (cm"^(K - km s-^)"^) = 
6.3 X 10" X aco(M0 pc~^(K - km s"^)"^) 



(2) 
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Hereafter, we refer to the CO-H2 conversion factor in terms of 
^ccEl though plot in terms of both Xco and aco • 

Despite potential variations in CO abundances, radiative trans- 
fer effects, and varying H2 gas fractions, a variety of inde- 
pendent measurements of H2 gas mass in GMCs have shown 
that the CO-H2 conversion factor in Galactic clouds is reason- 
ably constant, with Xco ~ 2 - 4 x 10^" cm~^/K-km s"^ 
(aco « 3 - 6 Mq pc"^ (K-km s"^)"^). Methods for obtain- 
ing independent measurements of H2 gas mass include (i) as- 
suming the GMCs are in virial equilibrium, and utilising the 
CO line widt h to derive the H2 mass within th e CO emitting 
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1986: Bertschetal. 1993: Strong & Mattox ' 19961 : iHunter et al] 
1997; Abdo et al.. . 2010b : ,Delahave et al. 201 1|) . Beyond this, ob- 
servations of GMCs in the Local Group suggest that a simi- 
lar CO-H2 conversion factor may appl y for some c louds out- 
side of our own Galaxy jRo solowsky ^tal] l2003l : iBlitz et all 
I2007I : Id onovan Meyer et al][201 If) . Numerical models of molec- 
ular clouds on both resolved and galaxy-wide scales have indi- 
cated that a relatively constant CO-H2 conversion factor in the 
Galaxy and nearby galaxies may naturally arise from GMCs that 
have a limited ra nge in surfac e densities, m etallicities and ve- 
locity dispersions jGioveT & Mac Low! [20 111: Ishetty et al]|2011bl : 
iNaravanan et alj201 lbl : lFeldmann et al J2OI ih . 

In recent years, a number of observational studies have pro- 
vided evidence for at least two physical regimes where the CO- 
H2 conversion factor departs from the "standard" Milky Way 
value. The first is in high-surface density environments. In- 
terferometric observations of present-epo ch galaxy mergers b y 
Scoville et al. ( 1991); Downes et al. ( 1993 l);ISolomon et alH1997l); 
bownes & Solomon (2003); Hinz & Rieke ('2006') : 'M eier etal.l 
(2010|) 

and Downes & Solomon ( 1998) showed that using a Milky 



Way Xco would cause the inferred H2 gas mass to exceed the dy- 
namical mass of the CO-emitting region for some galaxies. This 
implies that the CO-H2 conversion factor should be lower than 
the Galactic mean in high-surface density environments. More re- 
cent observat i ons o f z ~ 2 Submillimetre Galaxies (SMGs) by 
iTacconi et alj ( I2OO8I) suggested a similar result for high-redshift 
starbursts. Similarly, observations of GMCs toward the Galactic 
Centre indicate t hat Xco may be lower in this high-surface den- 
sity environment foka et alj 19981) . 

Second, in low-metallicity environments at both low and 



^ In the literature, Xqq is sometimes referred to as the X-factor, and we 
shall use the two interchangeably. 



.Leroy et al. 201 1 ; Genzel et al, 201 la), though t here is some debate 
ov er this (see summar ies in iBlitz et alj|200l and the Appendix 
of ITacconi et"ai] ( |2008|) ). Observations have suggest ed that the X- 
factor may scale a s X oc (O/H)"'' where b = 1—2.7 ( I Arimoto et alj 
Il996l ; llsraellll997l) . 

The fact that these two effects drive Xco in opposite direc- 
tions complicates the interpretation of CO detections from high- 
redshift systems where galaxies display a large range in metallici- 
ties (e.g. Shapley et al. 2004; Genzel et al. 2011a; Shapley 20 ij) 
and gas surface densities (e.g. iBoth well et alj 120 101 ; IDaddi et^ 
l2010al : lGenzel et alJlioiol : INaravanan et alj|2011al) . Further mud- 
dying the interpretation of high-z molecular line emission is 
the fact that there are not always clear analogs of high-redshift 
galaxies in the present-day Universe. For example, relatively 
unperturbed discs at 2: ~ 2 oftentimes have surface densities, 
star formation rates , and velocity dispersion s comparable to lo- 
cal galaxy mergers (jPaddi et alj|200^ , l2010al : iKrumholz & Dekell 
'2010'; Genzel et al.'2011b'), though (sometimes) lower metallicities 
(Cresci et al. 2010). Similarly, even the most heavily star-forming 
galaxies at 2 ~ 2, SMGs, at times show dynamically cold molec- 
ular discs even when they are potential l y the result of mergers 
frjarayanan et al. 2009l: ICarilli et al.ll2010l : lNaravanan et ai]|2010bl: 
lEngel et al. 2010). Converting CO line intensity to H2 gas masses 
is a multi-faceted problem that involves understanding how galactic 
environment sets the X -factor. 

Over the last two decades, models of GMC evolution have 
made substantive headway in elucidating the variation of Xco 
with physical properties of molecular clouds. The earliest GMC 
models u tilised ID radiative transfer with sp herical cloud mod- 
els (e.g iKutner & Leun3 1 19851 : IWalll 120071) . Photodissociation 
region (PDR) models furthered these studies by including the 
forma tion and destruction pathways of CO (Belletal. 200^ 
I2OO7I ; j Meiierink et al. I I2OO7I) . More recently, magnetohydrody- 
namic model s of GM C e volution with time-dependen t chemistry by 
iGlover et all ( I2OI0I) and lClover & Mac Lowl (l201lh coupled with 
radiative transfer calculations I ShettvetalJl2011bllj) have investi- 
gated Xco on the scales of individual GMCs, and its dependence 
on the physical environment. 

Compared to models of isolated GMCs, there are relatively 
few simulations explo ring the effect of the larg er galactic environ- 
ment on the X-factor. lMaloney & BlacM ( |l9 88') presented some of 
the earliest models which explored the effects of changing individ- 
ual physical para meters in isolation on th e CO-H2 conversion fac- 
tor. V ery recently , 'Feldmann et alJ ( I2OI ll) have tied the GMC mod- 
els of iGlover & Mac Low (2011|) to cosmological simulations of 
galaxy evolution to investigate Xco on galaxy-wide scales in rela- 
tively quiescent disc galaxies. Buildin g on these models, as well a s 
what has been learne d in the studies of IGlover & Mac Low! (1201 lb. 
IShettv etalJ ( l2011bl) and lWolfire et al.( ( I2010I) . in lNaravanan et alj 
( l2011bl) . we combined dust and molecular line radiative transfer 
calculations with hydrodynamic simulations of galaxies in evolu- 
tion in order to develop a model that aims to capture the CO line 
emission from GMCs on galaxy-wide scales. 

In this paper, util ising the methodology we developed in 
iNarayanan et alj jlQX Ibl) , we investigate the effect of galactic envi- 
ronment in setting the CO-H2 conversion factor in galaxies at low 
and high-redshifts. Our paper is organis ed as follows. In § [2| we 
summarise the methodology developed in lNaravanan et al] ( 201 Ibl) 
and employed here. We note that while § |2| is a shorter summary. 
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a more complete description is presented in Appendix [A] In § [3] 
we investigate the role of galactic environment on Xco, focusing 
on isolated disc galaxies (i) l3.1b . galaxies at low metallicity (ij |3.2b . 
high-surface density (§[33J, and high redshift (i; l3.4b . Building on 
these results, in § |4l we develop a functional form for calculating 
Xco from observations of galaxies and, as an application, utilise 
these to interpret Kennicutt-Schmidt star formation rate relations at 
low and high-redshift. In §[5] we discuss our results in the context 
of other theoretical models, and in § (6] we summarise our main 
results. 



2 SUMMARY OF SIMULATION METHODOLOGY 

Our main goal is to simulate the impact of galactic environment on 
the H2 content and CO emission from galaxies. This involves sim- 
ulating the evolution of galaxies, the physical state of the molec- 
ular ISM, and the radiative transfer of CO lines through GMCs 
and through galaxies. Because much of the methodology has been 
described in a previous paper by us ^Narayanan et alj 201 Ibh . we 
briefly summarise our approach here, and defer the quantitative de- 
tails to the Appendix. 

We first require model galaxies to analyse. We simulate the 
hydrodynamic evolution of disc galaxies in isolation and galaxy 
mergers over a range of galaxy masses, merger mass ratios, and 
redshifts utilising a modified version of the publicly available code, 
GADGET-2. These simulations provide information regarding the 
kinematic structure, mass and metal distribution of the ISM, as well 
as the stellar populations. Table AI summarises the model proper- 
ties. 

The remainder of our calculations occur in post-processing. 
We smooth the SPH results onto an adaptive mesh. We require 
knowledge of the physical and chemical state of the molecular 
clouds in our model galaxies. We assume the cold H2 gas is bound 
in sp herical GMCs with H2 f racti ons calculated fo l lowing the mod- 
els of lKrumholz et al.l ( l2008h and lKrumholz et al.l ( 12009 j) . Carbon 
is assumed to have a uniform abundance within these clouds of 
1.5 X 10"'' X Z' , where Z' is the metallicity with respect to solar. 
The fracti on of carbon lo cked up in CO is determined following the 
models of IWolfire et"ai] (|2010), and have an explicit dependence 
on the metallicity. GMCs that have a surface density greater than 
100 Mq pc~^ are considered resolved. GMCs that are not resolved 
(typically low mass GMCs in large cells in the adaptive mesh) have 
a floor surface density of the aforementioned value imposed, con- 
sistent with observations of Loca l Group GMCs (e.g. Bolatto et al, 
l2008l : lFukui & Kawamurj|20IOl) . Unresolved GMCs have velocity 
dispersions equal to the virial velocity of the GMC, whereas re- 
solved GMCs have velocity dispersions determined directly from 
the simulations. 

The temperature of the GMCs is determined via a balance of 
the various heating processes on the gas (here, photoelectric effect, 
cosmic rays and energy exchange with dust), and line cooling. The 
dust temperature is calculated via the publi cly ava ilable dust radia- 
tive transfer code, StJNRISE jjonsson et al. The cosmic ray 
heating rate is assumed take on the mean Galactic value (except 
as described in Appendix ICt. Utilising this model, the temperature 
ranges typically from ~ 10 K in quiesce nt GM Cs to > 100 K in 
the centres of starbursts l lNarav anan et al. 201 Ibh . 

Once the physical and chemical state of the ISM is known, 
we are prepared to model the CO line emission via radia- 
tive transfer calculations. The emergent CO line emission from 
the GMCs is calculated via an escape probability formalism 



( iKrumholz & ThompsonI |2007|) . This radiation is then followed 
through the galaxy in order t o account for radiative transfe r pro- 
cesses on galaxy-wide scales (Naravanan e t alj2006l . l201Ibl) . With 
these calculations, we know the thermal and chemical state of the 
molecular ISM in our model galaxies, the synthetic broadband 
SEDs, and the modeled CO emission line spectra. At this point, 
we are in a position to understand variations in the CO-H2 conver- 
sion factor. We remind the reader that further details regarding the 
implementation of these models can be found in the Appendix. 

We note that, in order to alleviate confusion between simula- 
tion points and observational data on plots, we will employ a sys- 
tem throughout this work in which filled symbols exclusively refer 
to observational data, and open symbols refer to simulation results. 



3 THE EFFECT OF OF GALACTIC ENVIRONMENT ON 

Our general goal is to understand how Xco depends on the phys- 
ical environment in galaxies, and how it may vary with observable 
properties of galaxies. In order to do this, we must first develop in- 
tuition as to how various physical parameters affect Xco- In this 
section, we examine how Xco varies from the Galactic mean in 
low-metallicity environments, high-surface density environments, 
and at high-redshift. 

Quantitatively, we define the mean Xco from a galaxy as: 



(Xco) = 



/ Sh2 dA 
/ M^co dA 



(4) 



which is equivalent to a luminosity-weighted Xco over all GMCs 
in the galaxy. 



3.1 Review of Xco in ^=0 Quiescent Disc Galaxies 

We begin by considering Xco in disc galaxies at z=Q with metal- 
licities around solar {Z' « 1). These galaxies have mean X-factors 
comparable to the Galactic mean and will serve as the "control" 
sample from which we will discuss variations in physical parame- 
ters. 

Recalling § |2] and referring to the Appendix, when clouds are 
not resolved, we impose a floor surface density of 100 Mq pc~^. 
In these GMCs, the velocity dispersion is the virial velocity of the 
cloud. In our model 2=0 disc galaxies, gas compressions within the 
galaxy are unable to cause significant deviations from these subres- 
olution values of Xh2 ~ lO'^^cm^^ and a ~ 1 — 5 km s"^. In this 
regime, the temperatures of the GMCs typically fall to ~ 8 — 10 K. 
This is the usual temperature where cosmic rays dominate gas heat- 
ing in our model; the densities are not sufficiently high for any of 
the heating processes to increase the gas temperatures drastically. 
At ~solar metallicities, a sufficient column of dust can easily build 
up to protect the CO from photodissocation, and most of the car- 
bon in molecular clouds is in the form of CO. With these modest 
conditions in the clouds, and little variat ion throughout the galaxy 
(see Figure 2 of iNaravanan et alj|2011bl) , the modeled Xco tends 
to be ~ a few xlO^" cm~^/K-kms~^ (i.e. similar to the Milky 
Way mean), with the only notable exception being GMCs toward 
the galactic centre (Naravanan et al. 201 Ib)^ 

As pointed out by INaravanan et al.l ( l20Ilbl) , while various 
subresolution techniques are folded into our model disc galaxies, 
that we see relatively little variation in the GMC properties in our 
model 2=0 discs is a statement that the galactic environment is 
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Figure 1. Xqq versus mass-weighted mean metallicity (in units of solar) 
for all 2=0 model galaxies witli Sh2 ~ 100 Mq pc~^. Tlie contours rep- 
resent tlie number of snapshots in a given Xqq-Z' bin, with the numbers 
increasing with increasing lightness of the contour. The dashed line outer 
contour encompasses all model galaxies, regardless of thek gas surfac e 
density. Overlaid are o bservational data poi nts from Bol atto et alj j2008l) ; 
iLerov et al.l )201 ll) and lGenzel eTal] jloT la) . The solid line shows our best 
fit to the simulations and is expressed in Equation[8]and described in §|4] 



not extreme enough to cause significant deviations from the de- 
fault surface densities and velocity dispersions. The temperatures 
and velocity dispersions are allowed to vary freely with galactic 
environment, an d the surface den s ities have a floor v alue similar 
to actual GMC s telitz et al1l2007l : iBolatto et ahllioosh . As shown 
bv lShettvetai] ( l2011bl) . GMCs with physical parameters compa- 
rable to those observed in Galactic GMCs exhibit Xco values 
close to the Galactic mean, Xco ~ 2 - 4 x 10^"cm~^/K-km s"^ 



(aco ~ 3 — 6 Mq pc~^ (K-km s~^)~^). As we will show, in 
mergers in the present Universe, low metallicity galaxies, and in 
some cases, high-redshift discs, the physical properties of GMCs 
vary sufficiently that this is no longer the case. 

3.2 The Effects of Metallicity on Xco 

The metallicity of the gas can have a strong effect on the CO-H2 
conversion factor. In met al-poor gas, it is possible to have "CO- 
dark" molecular clouds jPapadopoulos et al] I2OO2I : IWolfire et al.l 
l201Ct) . In these regions, H2 can self-shield to protect itself from 
photodissociati ng UV radiation, whereas CO cannot and requires 
dust to survive jStemberg & Dalgamoll99^ : lHollenbach & Tielenj 
[1999.) . In these cases, we expect a larger fr action of CO-dark clouds 
and a rise in the CO-H2 conversion factor ( Malonev & Blackl 19881: 



IWolfireetal.ll2010l ; IShettv et al.ll2011 j ; iFeldmann et al.ll201lh . 

In Figure[T] we plot the emission-weighted mean Xco for our 
model 2=0 isolated disc galaxies as a function of metallicity with 
filled contours. To control for surface density effects, we plot mod- 
els with Eh2 ~ 100 Mq ~^ . The contours indicate the number of 
snapshots in a given Xco~^' bin. The outer dashed contour de- 
notes the contour that encompasses all of our models, at any gas 
surface density. 

The lowest metallicity model galaxies all have relatively high 
mean X-factors compared to the ~solar metallicity galaxies. The 
lowest metallicity models can have mean X-factors approximately 



Figure 2. Mean Xco versus mass-weighted mean surface density for 2=0 
model mergers when they are undergoing a starburst (e.g. the snapshots with 
the peak SFRs). The plotting symbols are shown in the legend, though a 
given symbol may be replaced by a square if the metallicity is lower than 0.5 
. The purple filled circles represent the compiled data from Tacconi et 5] 
1 2008). The black solid hne is the best fit for all simulation snapshots from 
Equation llOl and is discussed in §|4] 



an order of magnitude greater than the Galactic mean. As the galax- 
ies evolve and the ISM becomes enriched with metals, the carbon is 
almost exclusively in the form of CO, and the X-factor decreases. 

The predicted Xco-^' relation in Figure [T] matches well 
with recent observational data. Overlaid, w e plot recent result s 
from resolved regi ons i n nearby galax i es fromlSolatto et aL I( l2008h . 
iLerov et alj ^20 11*) and 'Ge nzel etalj ( 1201 lal) . Note that Sh2 (and 
thus Xco) in the observations depends on an estimate of the gas 
mass, which is obtained with various methods for these studiefl 
Both the observations and models show an upward trend in Xco 
with decreasing metallicity. While we defer a discussion of fitting 
Xco in terms of Z' to § |4] we denote our best fit model (discussed 
in Equation[8} by the solid line in Figure[T] 

We should note that the contours in Figure[T]indicate the range 
of possible X-factors at a given metallicity as returned by our mod- 
els. They are not cosmological, so do not connote any particular 
probabilities. The simple fact that observed data points lie within 
the contours with a similar Xco-^' trend suggest a reasonable 
match between our models and observed galaxies. When examin- 
ing the outer contour that encompasses all of our models (and not 
just those at a given surface density as in the greyscale contours), 
it is clear that there is a significant dispersion in Xco on either 
side of the best-fit line in Figure[T] Moreover, there is a significant 
dispersion in the observed data at a given metallicity (especially 
noticeable near Z' « 1). This suggests that there is a second pa- 
rameter controlling the X-factor in galaxies, aside from metallicity. 



2 iBolatto et al] )2008l) assume clouds are viri ahsed; if instead clo uds are 
marginally bound, then Xco would decrease. iLerov et al.l i201 ll) use in- 
frared emission to derive dust masses, and HI observations to derive a dust- 
to-gas ratio . If this dust to gas r atio does not map to the H2 gas, the derived 
Xco will change accordingly. I Genzel et al.N2010l) assume a scaling rela- 
tion SsFR Sgag. If a steeper relation were adopted (c.f. 14. 21 . then 
the values of Xco corresponding to high values of SgpR would decrease 
somewhat. 
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As we will show in the next section, this is the dynamical and ther- 
mal state of the molecular ISM. 



3.3 High Surface Density Galaxies 

3.3.1 Large Temperatures and Velocity Dispersions 

On average, in regions of high surface density, the emission- 
weighted mean Xco in a galaxy is lower than the Galactic mean 
value of ~ 2 — 4 x 10^° cm~^/K-km s~^. This has been shown 
observationally bv Tacco ni et alj ( l2008h . and seen in the models of 
[Naravanan et al. (2011b). 

By itself, an increase in surface density does not cause Xco 
to decrease. Rather, the opposite is true: per Equation[T] at high sur- 
face densities, Xco would increase if Wco were fixed. However, 
in high-EH2 regions, the velocity-integrated line intensity, Wco in 
fact increases even more rapidly than Eh2, causing a net decrease 
in Xco - To see why this is, consider the physical processes that are 
typically associated with increased surface densities. 

First, regions of high s urface den s ity are associated with 
higher star formation rate s iKennicutt 1998a; Krum holz et al.l 
l2009bl : lOstriker et alllioiol : [Ostriker & Shettv 201 1) . While a va- 
riety of theories exist as to why th is is the case (see a review of 
some of these models in Pfai]|2010h . in our simulations this is due 
to the fact that we explicitly tie our star formation rates to the 
volumetric gas density on small scales. With high star formation 
rates come hotter dust temperatures as the UV radiation heats the 
nearby dust grains. When the gas densities are > lO^cm^'^, the 
dust and gas exchange energy efficiently, and the gas te mperature 
appro aches the dust temperature (Goldsmith 2001 ; Juvela & Ysard 
l201lh . Hence, in regions of high surface density, the increased 
dust temperature driven by the higher star formation rates causes 
an increase in the gas kinetic temperature. Because the CO (J=l- 
0) line is thermalised at relatively low densities, the brightness 
temperature of the line is comparable to the kinetic temperature 
of the gas, and is thus increased in regions of high surface den- 
sity. For low redshift galaxies, the easiest way to increase the sur- 
face density tends to be through merging activity and the associ- 
ated t idal torques which d rive gaseous in flows into the nuclear re- 
gion jBames & Hemquilt 1991. 1996; Mihos & Hemauisllll994bl . 
1 19961) . though higher-redshift discs can have relatively large sur- 
face densities si mply due to gra vitational instabilities in ext remely 
gas-rich clumps jSpringel et al . 2005; Bo umaud et alj2oTol) . 

Second, in the simulations, high surface-densities in gas are 
typically accompanied by high velocity dispersions. This gener- 
ally means some level of merging activity in low-redshift galax- 
ies, and either merging activity or unstable gas clumps in high- 
redshift discs. In major mergers, the emission-weighted velocity 
dispersions can be as high as --^ 50 — 100 km s~^. 

The increased velocity dispersion and kinetic temperature 
contribute roughly equally to increasing the velocity-integrated line 
intensity. During a merger, the temperature and velocity dispersion 
individually increase enough to offset the increase in the gas surface 
density, and in combination tend to drive the emission-weighted 
mean Xco for a galaxy below the Milky- Way average. In regions 
of high surface density, the emission-weighted mean Xco tends 
to decrease on average. In 13.41 we detail specific numbers for a 
sample merger that serves as an example for this effect. 

We can see this more explicitly in Figure|2l where we plot the 
mean Xco versus mean surface density for the starburst snapshots 
in our low-redshift sample of galaxies. The starburst snapshots are 
defined as the snapshots where the SFR peaks for a given model 



and are categorised by the merger mass ratio (1:1, 1:3 and 1:10). 
The mean surface density is defined as the mass-weighted surface 
density over all GMCs, i, in the galaxy: 

. , _ E. SH2.i X MH2.i 

Henceforth, when we refer to the surface density of the galaxy, we 
refer to the mean surface density defined by Equation [5] We note 
that this is different from the commonly used definition Eh2 = 
Mu2/A, where A is the area within an observational aperture. We 
refrain from the latter definition as it is dependent on the choice of 
scale. (Eh2) can be thought of as the surface density at which most 
of the mass resides. Immediately, we see two trends in Figure[2l 

First, with increasing surface density, we see decreasing mean 
Xco due to the warm and high velocity dispersion gas associated 
with merging systems. Second, the most violent mergers tend to 
have the lowest mean Xco values, whereas lower mass ratio merg- 
ers (1:3) have less extreme conditions, and thus X-factors more 
comparable to the Galactic mean v alue. The purple circ les in Fig- 
ure[2]note observational points fror n l Tacconi et al]i2008l) ; the mod- 
els and data show broad agreemenjj. 

The open squares in Figure|2]are points from our models with 
mean metallicities less than Z' < 0.5. Careful examination of these 
points shows that some of them have rather large X-factors. As 
we discussed in § 13.21 lower metallicity galaxies have larger mass 
fractions of CO-dark clouds, and thus higher X-factors. 

In Figure [3] we demonstrate more explicitly the effect of 
metallicity on the Xco- surface density relationship. Because we 
find it useful in a forthcoming section (§ ^ to parameterise Xco 
in terms of the observable velocity-integrated CO intensity, H^co 
as an observable, rather than Sh2, we plot Xco against H^co in 
Figure[3] H^co is defined as the luminosity-weighted line intensity 
from all GMCs in the galaxy. We plot the Xco -Wco relationship 
for all snapshots of all 1:1 2=0 major mergers within two distinct 
metallicity ranges. The selected metallicity ranges are arbitrary, and 
are chosen simply to highlight the influence of metallicity. As is 
clear, the lowest metallicity points in Figure[3]have the highest X- 
f actors, and the reverse is true for the highest metallicity points. 
For each metallicity bin, the trend is such that increasing Wco (or 
Eh2) correlates with decreasing Xco. though the normalisation 
varies with metallicity. This informs our fitting formula in §|4l 

Returning to Figure (2] we make a final point that there is a 
large dispersion in Xco for the merger models. Both the 1:3 mod- 
els and 1:1 mergers show Xco values ranging from above the 
Milky- Way mean to an order of magnitude below it during their 
peak starburst. When sampling the entire library of merger orbits 
for a given merger mass ratio, a wide range in outcomes is appar- 
ent. We see a diverse set of velocity dispersions in the gas, as well 
as star formation rates, owing to differing efficiencies at which an- 
gular momentum is removed from the gas. Some models undergo 
rather vigorous starbursts (approaching ~ 500 MQyr~^), whereas 
others hardly sustain a noticeable starburst upon final coalescence. 
Galaxies that undergo their peak starburst only on first passage can 
have rather different metallicities in their ISM than mergers that go 

We caution that our definition of (Sh2 ) is a mass- weighted surface den- 
sity, and is different from the surface density defined by the Tacconi et al. 
data, Mh2/j4-. In the limit of a large volume filling factor, these values will 
approach one another. The observed data closer to Sn^oi = 100 Mq pc~^ 
likely represents galaxies with a clumpy ISM, and the modeled and ob- 
served H2 gas surface densities may differ in this regime. 
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Figure 3. Xco versu.s CO line intensity, Wqo for 2=0 1:1 galaxy mergers 
in two distinct metallicity bins. In this case, Wco is a suiTogate for Sh2 
that is an actual observable (and relateable via Xco)- At a given metalUcity, 
Xco decreases with CO intensity due to the larger number of CO-dark 
GMCs in these galaxies. The solid line shows the best fit model (c.f. § |4) 
for a metalHcity of Z' = 1. 

through a vigorous star formation period during first passage and 
inspiral before experiencing a starburst contemporaneous with fi- 
nal coalescence. 

In the solid line of Figure [4] we plot the distribution of 
emission-weighted mean Xco values for each of our 1:1 z=Q 
merger models during the peak of their star formation rate. For 
comparison, we plot the distribution of ^-factors for our z=Q discs 
in the green dashed line (only plotting galaxies with ~solar metal- 
licity), as well as the distribution of Xco for our z=2 discs (which 
we discuss in more detail in ii l3.4t . As we see, there is no "merger" 
value for Xco - It is possible to have Xco values in starbursting 
mergers comparable to the Milky Way's. The fact that starbursting 
mergers, on average, have lower X-factors than the Galactic mean 
is likely due to a selection effect. We return to this point in more 
(quantitative) detail in § |4] 

3.4 Xco in High-Redshift Galaxies 

3. 4. 1 Basic Results 

Now that we have developed intuition regarding the variation of 
Xco with metallicity and surface density, we are in a position to 
understand how galaxies at high redshift may behave with respect 
to the CO-H2 conversion factor. 

Mergers at high-z are some of the most luminous, rapidly 
star-forming galaxies in the Universe. As an example, many 2: ~ 2 
Submillimetre-se lected G alaxies (SMG) form stars at > 1000 
Mq /yr l lNaravanan et alj 120109) . However, despite the ~order 
of magnitude greater star formation rate in these galaxies com- 
pared to local mergers, mergers at high- 2: are similar to their low- 
redshift counterparts in terms of their typical X-factors. While the 
mean gas surface densities in e.g. z ~ 2 Submilli metre Galaxies 
(SMGs) are larger than low- redshift mergers (e.g. iTacconi et al.l 
l2008l : [Narayanan et alj|2010bh . both the dust temperatures and gas 
velocity dispersions also rise commensurately. 

As a specific example, we focus on a 1:1 major merger 



Figure 4. Distribution of Xco values during the peak of the star formation 
rate (the "burst" snapshot) for all 1:1 and 1:3 local galaxy mergers (black 
solid), all z=0 1:1 and 1:3 mergers (i.e. not just the burst snapshots; black 
dash-dot), all high-z disc models with Z' 1 (red dotted), and low-z disc 
models with Z' ^ 1 (green dashed). See text for details. 

at z ~ 2 (This is model "z3b5e"; please refer to the Table in 
the Appendix for the initial conditions of this model merger). 
Model z3b5e undergoes a luminous burst of ~ 1500 MQyr"^, 
and may be selected a s a submiUimetre gala xy when it merges 
jNaravanan etall l2009l : iHavward et J] I2OIII) . During the burst, 
this simulation reaches a mass-weighted mean surface density of 
molecular gas of ~ 10* Mq pc~^. At the same time, the mass 
weighted kinetic(dust) temperature is ~ 150(1 60) kQ, and the mass- 
weighted velocity dispersion in the GMCs is ~ 140 km s^^. Do- 
ing a simple scaling results in a mean Xco of ~ 5 x 10^'"* 
cm^^/K-km s^^ . Of course the real mass-weighted value may vary 
from this owing to both radiative transfer as well as the fact that this 
simple scaling is not a true averaging. Observational estimates of 
the X-factor in SMGs suggest that they are similar to the lower 
values in the range of local ULlRGs ( ITacconi et alj2008h . 

Our simulated disc galaxies at high-redshift show a range of 
Xco values, ranging from comparable to the Galactic mean to val- 
ues 2-5 times lower. The reason massive discs at high-redshift may 
have lower X-factors than the Galactic mean can be understood 
in the following way. In contrast to present-epoch galaxies, galax- 
ies at higher-redshifts (z>_l)_ata fixed stellar mass are denser 
and more gas-rich (e.g . |Erb et al. '2006'; Forster Sc hreiber et al.l 
,200ft ; iDaddi et alj|2010bl ; ITacconi et al. 201Q ). Both simulations 
and observations suggest that galaxies around redshifts z ~ 1-2 
may have baryonic gas fractions of order 20-60% ( iDave et alj2"oiol ; 
Daddi et al. 2010b; Tacconi et al. 2010). A primary consequence of 
this is that discs at higher-redshifts may be heavily star-forming, 
with star formation rate s or order ^ 10^ Mf^yr^^, comparable 
to local galaxy mergers jPaddi et al]l2007l ; lForster Schreiber et al.l 
l 2009| ; lNaravanan et alJ2010 a^ In fact, simulations suggest that disc 

S ;alaxies at z ^ 2 lik ely dominate the infrared luminosity function 
Hopkins et alj201(t) . 

* Note that owing to radiative transfer effects, this dust temperature is not 
necessarily what would be derived simply by identifying the location of the 
peak of the SED. 
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In the absence of rather extreme stellar feedback, very gas- 
rich discs at high-redshift can be unstable to fragmentation , 
and form massive ^ kpc-scale clumps (e.g. Springel et al .1 1200^ : 
ICeverinoetai]|201(3[) . These clumps can have relatively high ve- 
locity dispersions (~ lO'^ km s~^) and warm gas temperatures 
owing to high volume tric densities and high star formation rates 
teoumaud et alj2010h . 

These effects are the strongest in the most massive discs. Our 
most massive 2=2 model disc galaxy has a total baryonic mass of 
~ 5 X 10^^ MQyr~^, and has typicaljX-factors ranging anywhere 
between a factor of five below the Galactic mean to the Galactic 
mean value. The lower mass z=2 isolated disc models (with bary- 
onic masses of Mbar = 1 X 10" and 3.5 x lO"' Mq) typically 
have X-factors comparable to the Milky Way mean. Returning to 
Figure [4] we examine the red dotted line that represents z ~ 2 disc 
models. Because the idea of a 'starburst' snapshot is less meaning- 
ful for the evolution of a disc galaxy, we plot the X-factor for every 
snapshot for our model discs with metallicities around solar. We see 
a large spread in mean X-factors. 



3.4.2 Do Mergers and Discs Have Inherently Dijferent 
X -factors? 

In light of the fact that high-z discs have, at times, star formation 
rates comparable to local galaxy mergers, a pertinent question is 
whether there is an intrinsic difference in the X-factor between 
high-z discs and galaxy mergers. Another way of saying this is, 
for a given set of physical conditions, are the X-factors from merg- 
ers lower than the X-factors from high-z gas-rich, gravitationally 
unstable discs? A cursory examination of Figure |4] indicates that 
mergers (the black solid line) have systematically lower X-factors 
than discs (the blue and red dashed and dotted lines). Indeed, in 
the local Universe, it is obser ved that mergers ha ve, on average, 
lower X-factors than discs (e.g. lTaccoiri et alj2008h . However, this 
is likely due to a selection bias. We remind the reader that the black 
solid line in Figure|4]represents starbursting mergers. These merg- 
ers are caught when their gas is extremely warm and with large 
velocity dispersion. When comparing mergers and discs with com- 
parable physical conditions, the observed Xco values are in fact 
quite similar. It is the physical conditions in a galaxy that deter- 
mine the X-factor, not the global morphology. 

To demonstrate this, we perform three tests. First, we compute 
the distribution of Xco values for all \ :\ and 1 :3 merger snapshots 
(at ~solar metallicity), and indicate this with the black dot-dashed 
line in Figure|4] As we see, the distribution of Xco values is broad, 
but with substantially less power in the low X-factor regime than 
the distribution that denotes only starbursting mergers (black solid 
line). This highlights that mergers which are selected during a par- 
ticularly active phase are more likely to have low X-factors, due 
to their warm and high-u gas. When controlling for this effect by 
picking galaxies with similar CO intensity (Wco) and metallicity, 
mergers and discs have the same Xco on average. 

^ Because our simulations are not cosmological, there is no accretion of 
intergalactic gas. As a result, the metallicities in our model galaxies only 
rise with time. Because the jf -factor is dependent on metallicity (§ 13. 2K 
we have to make a choice as to which snapshot/metallicity to consider as 
a 'typical' galaxy. We assume any snapshot above Z' > 0.5 is "typical" 
based on the steady-state metallicities found for galaxies of baryonic mass 
compai'able to those in our sample from cosmological modeling (Figure 2 of 
[pave et al. 2010, though see Keres et al. (201 1) ; Vogelsberger et al.. 1 2011.) 
and fSiiacki et alj <201ll) ). 
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Figure 5. Comparison of the X-factor between low-z mergers (1:1 and 1:3) 
and high-z star forming discs. The histogram denotes ratio of X-factor from 
mergers versus high-z discs between snapshots with a similar metallicity 
and CO intensity. The sharp peak near unity implies that galaxies with sim- 
ilar physical conditions have siinilar X-factors, independent of large-scale 
morphology. 



To show this, in Figure|5] we perform our second test in which 
we examine the X-factors from all the 1:1 and 1:3 mergers (at low 
z) and compare them to the Xco from high-2 discs with the sam^fl 
metallicity and CO intensity {Wco)- We could equivalently per- 
form this analysis in terms of Eh2, though as we will show in §|4] 
parameterising in terms of Wco is desirable with regards to ob- 
servations. There is a strong peak at Xco ratios near unity, with 
some spread. The median value in the distribution is ~ 0.8, and 
the mean is ~ 1.1. The implication from Figure [5] is that galaxies 
with similar physical conditions (here, Z' and Wco) have similar 
X-factors, regardless of whether they are discs or mergers. The fact 
that mergers, on average, have lower X-factors than discs in the lo- 
cal Universe likely derives from the fact that they are selected as 
starbursts, which have preferentially higher temperatures and ve- 
locity dispersions in the gas. 

Third, in Figure|6] we examine the relationship between Xco 
and Wco for the same galaxies plottecQin Figures |4] and [5] These 
are all galaxies with metallicities around solar. The principal result 
from Figure [6] is that galaxies within a relatively limited metallic- 
ity and Wco (or surface-density) range have similar X-factors, re- 
gardless of the type of merger it is. Mergers and discs have similar 
Xco values when they have similar physical conditions, and are 
not inherently different based on their global morphology. In addi- 
tion, Figure[6l like Figures|2]and[3] shows a systematic decrease of 
Xco with increasing Wco (and Smoi)- 



^ "The same" here means that the values of Z' and Wco ^re within 10% 
of one another. 

^ To reduce clutter in the plot, we randomly draw 10% of the galaxies 
within each merger ratio bin to plot. 
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4 APPLICATION TO OBSERVATIONS 
4.1 Deriving Xco from Observations 

As we have seen from the previous sections, it is clear that there is a 
continuum of Xco values that vary with galactic environment. The 
dominant drivers of the X-factor in our simulations are the metal- 
licity of the star-forming gas, and the thermal and dynamical state 
of the GMCs. Informed by this, we are motivated to parameterise 
Xco as a function of observable properties of galaxies. 

Metallicity is a crucial ingredient to any parameterisation. 
At subsolar metallicities, we see the rapid growth of CO-dark 
GMCs. This has been note d both in observations (e.g. lLerov et al.l 
I2OIII: lGenzeletalJ|2011ah. as well as other numerical models 
iaf;(li ' ■ 



dShettv et alj20 11al.lKru mholz et alj2011bl : lFeldmann et alj20l"Th . 

As we saw in 8^ 13. 21 as well as in Figure[3] at a given galaxy surface 
density (or CO intensity), Xco increases with decreasing metallic- 
ity. 

Beyond this , as was shown in §[33] as well as in the models of 
[Narayanan et alj ( l2011bh . galaxy surface density is correlated with 
the thermal and dynamical state of the gas: at a given metallicity, 
higher surface density galaxies, on average, correspond to galaxies 
with a warm and high velocity dispersion molecular ISM, due to 
their higher SFRs. 

Informed by these results, we perform a 2D Levenberg- 
Marquardt fit (Markwardt 2009) on our model galaxies (consid- 
ering every snapshot of every model), fitting Xco as a function of 
mass-weighted mean metallicity and mass-weighted mean H2 sur- 
face density. We find that our simulation results are reasonably well 
fitted by a function of the form: 
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Figure 6. Comparison of Xco versus CO intensity (Wco) for low-2 
galaxy mergers and high-2: discs in an effort to investigate if mergers and 
discs inherently have different Xco properties. Included in this plot are 
all 1:1 and 1:3 mergers simulated at 2=0. Only snapshots with metallic- 
ities Z' > 0.7 are shown. To reduce clutter in the plot, we plot only a 
randomly drawn subsample (10%) of the snapshots from each mass ratio. 
The line shows the best fit from Equation [8] Evidently, galaxies that have 
similar physical conditions have similar X-factors, independent of galaxy 
morphology or evolutionary history. See text for details. 



Xco 



1.3 X 10^^ 
Z' X {Eh2)°-5 



(6) 



where Eh2 is in units of Mq pc^^ and Xco is in units of 
cm~'^/K-km s~^. Equation |6] provides a good fit to the model 
results above metallicities of ^' ~ 0.2. Turning again to Fig- 
ures [T] |2l [3] and (6) we highlight the solid lines which show how 
Equation [6] fit s both the simulation resu lts and observational data. 
We note that iOstriker & Shettvl ( l201lh obtained a similar result, 
aco oc by interpolating between empirical qco values 

{o-co = 3.2 for Snioi = 100 Mq pc~^ and qco ~ 1 for 
E^oi = 1000 Mq pc-2). 

Because Eh2 is not directly observable (hence the need for an 
X-factor), we re-cast Equation|6]in terms of the velocity-integrated 
CO line intensity. In order to parameterise Xco in a manner that 
is independent of the effects of varying beam-sizes or observa- 
tional sensitivity, we define the observable CO line intensity as the 
luminosity-weighted CO intensity over all GMCs, i: 



{Wco) 



J Wco dA 



T,Lco,i X Wco.i 
J2Lco,i 



(7) 



where {Wco) is in units of K-km s ^, and is the CO surface bright- 
ness of the galaxy. We then fit to obtain a relation between Xco, 
Z',and {Wco): 



Xno ~ 



6.75 X 10^° X (Wco)' 

2; 10. 65 



(8) 



where again (Wco) is CO line intensity measured in K-km s^^ , 
Xco is in cm^^/K-km s~^, and Z' is the metallicity divided by the 
solar metallicity. By converting Xco to aco, we similarly obtain: 



aco 



10.7 X (Wco)- 



(9) 



Where qco is in units of Mq pc~ (K-km s~^)~ . 

It is important to recognize that the power-law in Equation [8] 
cannot describe Xco indefinitely. At very low Wco, GMCs tend 
toward fixed properties in galaxies and the galactic environment 
plays a limited role in setting Xco- Considering this. Equation [8] 
formally becomes: 



Xco 



min [4,6.75 x (Wco)'"'^^] x 10^ 

^/0.65 



or, similarly: 



aco = 



[6.3, 10.7 X (Wco)- 



(10) 



(11) 



Equations [To] and [TT] can be used directly with observations 
of galaxies to infer an expected X-factor. One advantage of this 
formalism is that it captures the continuum of CO-H2 conversion 
factors, rather than utilising bimodal "Disc" and "ULIRG" values. 
Because we have chosen the physical quantities in our modeling 
based on mass or luminosity-weighted averages, they are defined 
without reference to a particular scale. Consequently, Equation [To] 
can be used on scales ranging from our resolution limit of ~ 70 pc 
to unresolved observations of galaxies. 

It is conceivable that alternative definitions of the observed 
mean CO intensity could be appropriate. One can imagine imple- 
menting an area-weighted intensity, i.e. Wco = ico/Area. This 
has the undesirable attribute of being dependent on a defined scale. 
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4.2 The Kennicutt-Schmidt Star Formation Relation in 
Galaxies from z=(i-2 

A natural application of our model for the CO-H2 conversion fac- 
tor is the Kennicutt-Schmidt star formation rate surface density- 
gas surface density relation in galaxies. Because the inferred H2 gas 
masses from observed galaxies are inherently dependent on conver- 
sions from CO line intensities, our understanding of the Kennicutt- 
Schmidt relation is fundamentally tied to the potential variation of 

Xco with the physical environment in galaxies^ 

Recent surve ys of both local galaxies (e.g. lKennicuttlll998al lbl 
iBigiel et al.l2008l and referenc es therein), as well as pioneering ef- 



forts at higher redshifts ( e .g. Bouche et alj 20071 ; iBothwell et alj 



l20ld ; iDaddi et alj l2010j ; lOenzel et alj l2010h have provided a 
wealth of data contributing to our knowledge of the star for- 
mation relation in both qu iesce nt disc galaxies an d starbursts. 
Work bv lDaddi et id] fcOlOiJ) and lOenzel et^ ilOld) demonstrate 
the sensitivity of these relations to the CO-H2 conversion factor: 
When applying the traditional bimodal conversion factor (Xco~ 
8 X 10^^ cm~^/K-km s"^ for ULIRGs and Xco~ 2 x 10^" 
cm~'^/K-km s~^ for discs) to the starburst galaxies and discs, re- 
spectively, a bimodal SFR relation becomes apparent when the data 
are plotted in terms of Esfr and Emoi. 

In Figures[7]--[9l we illustrate the effects of our model fit for 
Xco as a function of galaxy physical properties. In the left panel 
of Figure |7] we plot the star formation rate surface density for 
both loc al ga laxies and high- 2: ga laxies as compiled bv lDaddi et alj 
( l2010ah and lOenzel et aU ^2010*) against their CO line intensity, 
Woo. Although there is significant scatter, the Esfr, vs W^co rela- 
tion is unimodal. In the middle panel, we plot the SFR-Eh2 relation 
utilising the bimodal X-factors assumed in the literature (with the 
above "ULIRG" value for the inferred mergers [local ULIRGs and 
high-z SMGs], and the above "quiescent" value for low-z discs and 
high-z BzK galaxiesfl The circles are unresolved observations of 
disc galaxies at low and high-z, triangles and contours are resolved 
observations of local discs, and the squares are local ULIRGs and 
inferred mergers at high-redshift. When separate high and low val- 
ues of Xco are adopted for discs and mergers, a bimodal relation- 
ship between Esfr and Egas results, with power-law index ranging 
between unity and 1.5. 

In the right panel of Figure [7] we apply Equation [TO] to the 
observational data (assuming Z' = 1 for the galaxies). The scat- 
ter in the modified relation immediately tightens, and it becomes 
unimodal. To numerically quantify the reduction in scatter with the 
modified relation we examine the ratio of the maximum inferred 
Eh2 to the minimum for all points within the relatively tight SFR 
surface density range of [0.05,0.1] MQyr~^ kpc^^ for the centre 
and right panels of Figure [T] The scatter is reduced by approxi- 
mately a factor of 5. Within this Esfr range, no mergers are in the 
sample. Thus, the reduction in scatter is not due to simply using 
a unimodal Xco versus bimodal Xco. We note that in applying 
Equation[TO]to the observed data, we have to assume that the inten- 
sity within the reported area is uniform. If the emission is instead 
highly concentrated over very few pixels, then the application of 
Equation [To] may overestimate Xco • 

The reason for the transition from a bimodal to unimodal KS 
relation is clear. In the modified relation, similar to the traditional 



* In practice, for the high-z BzK galaxies, we utilise a n X-factor of Xco ~ 
2.3 X 10^° cm"2/K-km s~^ to remain consistent with lOaddi et aljj20103) . 
though the usage of this versus the more standard disc value makes little 
difference. 



Kennicutt-Schmidt plot that uses a bimodal X-factor, the lower lu- 
minosity discs have CO-H2 conversion factors comparable to the 
Galactic mean, and the most luminous discs have X-factors up to 
an order of magnitude lower. Very massive, gas-rich, unstable discs, 
as well as lower luminosity mergers have X-factors in between the 
two, however, and fill in the continuum. Utilising Equation [TO] a 
simple linear chi-square fit of the observed data on the right side of 
Figure [7] returns: 



logio(EsFR) = 1.95 X logio(E 



mol / 



4.9 



(12) 



where Esfr is measured in MQyr ^ kpc and Emoi is mea- 
sured in MfT) pc~^. Utilising an em pirical method to obtain Xco oc 
^ lOstriker & ShetT 



iingj 

3 j20Ilh previously sh owed that the ob- 
servational data compiled in lGenzel et al. yields a similar fit 



to Equationll2l 

We remind the reader of the assumptions that have gone into 
this fit: We have assumed that every galaxy has solar metallic- 
ity, and neglected any potential effects of differe ntial excitation in 
the CO as a function of infrared luminosity (e.g. [Narayanan et alj 
l20Ilal) . Nevertheless, the application of a variable X-factor on the 
SFR-gas surface density relation has interesting implications. 

First, the index of ~ 2 of Equationll2lis consistent with the an- 
alytic models and hydrodynamic simu lations of self-regulated star 
formation by lOstriker & Shettvl ( 1201 II) . This work suggests that 
in molecular regions where supernova-driven turbulence controls 
the SFR and gas dominates the vertical gravity, the SFR surface 
density should be proportional to the gas surface density squared: 
log(EsFR) ~ 2 X loK(Emoi) — 5.0 (adopting fiducial parameters 
in Equation 13 of lOstriker & Shettvl ( I2OIII) ). This is shown as the 
solid line in the right panel of Figure [7] and is very comparable to 
the b est fit relation. Second, co mparing Equation [T2| to Equation 
13 of lOstriker& Shettvl ( I2OI l|), the empirical results are consis- 
tent with a value of momentum injected/total stellar mass formed 
of /p X p./m, ~ 3000 km s~^ (Emo i/lOO Mq pc-^)'°-°^; the 
fidu cial value adopted irilOstriker & Shettv ( 201 1) is 3000 km s~\ 

iDaddi et alj i20 lOal) and lGenzel et all 12010) suggest an alter- 
native mechanism for reducing the scatter imposed by the utilisa- 
tion of a bimodal Xco in FigurejT] Specifically, these authors find 
that by dividing the molecular gas surface density by the galaxy's 
orbital time, the observed Kennicutt-Schmidt relation goes from 
bimodal to unimodal, suggesting that the galaxy's global proper- 
ties are related to the local, small-scale processes of star formation. 
That is to say, when using a bimodal Xco, the Esfr — Emoi re- 
lationship is bimodal, whereas the Esfr — Emoi/idyn relationship 
is unimodal, with some scatter. 

If one abandons the bimodal Xco approximation, and utilises 
our favoured model for Xco, the observed relationship between 
Esfr and Emoi/idyn remains unimodal, and in fact the scatter in 
the relation is reduced compared to what one obtains using a bi- 
modal Xco . To show this, in FigurejS] we show the analog to Fig- 
ure |7] but with the abscissa showing the surface density or Wco 
divided by the dynamical time. The dynamical times used are the 
same as those in Daddi et al. (2010a) and Genzel et al. (20IG), and 
are defined as the rotational time at either the galaxy's outer radius, 
or half-light radius, depending on the sample. In the left panel of 
FigurejS] we show the relationship between Esfr and W^co/idyn 
(i.e. pure observables); in the centre panel, we show the relationship 
between Esfr and Emoi /tdv n wh en assuming a bimo dal Xco (as 
is done in lOaddi et al. I ( l2010j) and lGenzel etal] ( I2OI0I) ), and in the 
right panel we show the same relationship, but with Emoi deter- 
mined using our be st-fit continuous Xc o, ra ther than the bimoda l 
Xco value used in iDaddi et alj ( l2010ij) and iGenzel eTai] ( l2010h . 
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Figure 7. Kennicutt-Schmidt star formation relation (SFR surface density versus H2 gas surface density) in observed galaxies. Circles and triangles are local 
discs or high-^ BzK galaxies, and squares are inferred mergers (local ULIRGs or high-2 SMGs). Colours denoting separate surveys are described below. Left: 
SFR surface density vs. velocity-integrated CO intensity, yielding a unimodal SFR relation. Centre: When applying an effectively bimodal Xqq (oco =4.5 
for local discs, 3.6 for high-z discs, an d 0.8 for mergers), t he resulting SFR relation is bimodal. The solid and dotted lines overplotted are the best fit tracks 
for each "mode" of star formation as in lDaddi et alj ilOlOah . Right: SFR relation when applying Equation |10| to the observational data, resulting in a unimodal 
SFR relation. The power-law index in the relation is approximately 2 (solid line). Symbol legend: We divide galaxies into 'disc-like' with filled circles, and 
'merger-like' with squares. This assumes that high-z Bz^f galaxies are all discs, high-z SMGs and low-z UL IRGs are all mergers. T he low-z disc observa tions 
(black filled circles and black triangles) come from Kennicutt et al. (2007); Wong & Blitz (2002); Crosth waite & Turned i2007l) ; ISchuster et al.l j2007l) and 
comprise both resolved and unresolved points. The resolved points from the survey of Bigiel et al. (2008) are denoted by the coloured contours. The local 
ULIRGs are compiled by, Kennicutt ( 1998b) and are denoted by black filled squares. The high-z discs come fro m Ge nzel et al. (2010), Daddi et al. (2010b;J 
and are represented by filled blue circ l es. Th e high-z SMGs are divided into the the sa mples o f Bothwell e t al ] {20\&l (purple). iBouche et alj klOOt) (green), 
and lGreveerZI i2005l) ; lTacconi et alj j2003 . l2b08i) ; iEngel et ah C201Q) as compiled bv lGenzel et al.. C2010,) (filled red squares). 



We find a best fit relation (using our model for Xco) of: 

loglo(EsFR) = 1.03 X loglo(E„,ol/tdyn) ^ 1-05 (13) 

where Ssfr and Emoi/^dyn are both in MQyr^^ kpc~^. 

Our best fit Esfr — Smoi/idyn relation has slope of approxi - 
mately unity, comparable to what is found by Genzel et al. (ioitj), 
and consistent with the best fit slope of Daddi et al.,, (2010a,) of 
~ 1.15. A principal difference between using our model Xco ver- 
sus the bimodal Xco in calculating the Esfr — Smoi / idyn relation 
is a reduction of scatter. When measuring the scatter near SFR sur- 
face density of 1 Mq yr~^ kpc~^, we find that using our model 



Xco versus the bimodal value reduces the scatter by a factor of 
~ 5. 

The fact that our model for the Esfr — Smoi/idyn relation 
is consistent wit h the observed one ( thou gh with reduced sca tter) 
is not surprising. lOaddi et"al] l l2010ah and lGenzel et al.l dlOlOl) as- 
sume a Milky Way-like Xco for their disc galaxies, and roughly 
a factor of 5 lower for their mergers. In our model, the assump- 
tion of a bimodal Xco for discs and mergers is correct on average. 
The mean value for Xco for high-z SMGs and low-z mergers is 
in fact lower than the mean value for local discs (c.f. § 14. 3t . How- 
ever, many galaxies lie in the overlap region. Some local ULIRGs 
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Figure 8. Similar to Figure [T] but with the abscissa showing Wqo or ^mol divided by the orbital time of the observed galaxy. Symbols are the same as in 
Figure|2] but we omit galaxies for which orbital times are not available. The best fit slope in the right panel is of order unity. See text for details. 



have X-factors comparable to the Galactic average, and some high- 
z discs have X-factors more similar to the canonical literature 
"merger value". By modeling the continuous nature of Xco and 
more properly treating these intermediate cases, a reduction in scat- 
ter is natural. 

Finally, we consider Ssfr as a function of Emoi/iff, where 
tjj is the free-fall time within the dense molecular star-forming 
clouds in the galaxy. Both observations and theory have suggested 
that the star formation efficiency per free-fall time eg in molec- 
ular gas is approximately constant ( Krumholz & McKee 'loos'; 
Krumholz & Tan 2007; Evans et al. 2009; Ostriker & Shetty 201 1 ; 
Krumholz et al.ll201 laf) . We infer tg from the ob servable proper- 



ties of the galaxy (using the approximations of iKrumholz et al.l 
l2011ah , together with either a bimodal Xco or our favoured con- 
tinuous Xco. Krumholz et al. show that using tfi rather than tdyn 
makes it possible to fit the unresolved extragalactic o bservations, 
the re solved observations of Local Group galaxies from lBigiel et al.l 
( I2008 ). and individual molecular clouds in the Milky Way on a sin- 
gle relation, as illustrated in the left panel of Figure |9l The right 
panel of the Figure shows that this remains true for our favoured 



Xco and that a star formation law Esfr = effEmoi/ifi with 
efl ~ 0.01 remains a good fit to the observational data. As with 
the Ssfr — Smoi/tdyn relation, using our continuous Xco actu- 
ally reduces the scatter, and for the same reason: our favoured Xco 
produces essentially the same result as the traditional bimodal Xco 
for galaxies at the extremes of the disc and merger sequences, but 
makes the behaviour of Xco continuous rather than discontinuous 
for the large number of galaxies in the overlap region. 

The results f r om Fi gure [9] are compatible with the model of 
lOstriker & Shettvl ( |20T J). For a disc in vertical hydrostatic equi- 
librium with gravity dominated by the gas, and vertical velocity 
dispersion, v^, 

Esfr = effEsFR/tff- = etf x 4GELi/(\/3«z) (14) 

l lOstriker& Shettvl |20 111 Equation 21). Comparing to the fit ob- 
tained in Equation [12] using our continuous Xco relation, we 
find that Efi/u, = 0.001(km s^^ (En^oi/lOO Mq pc-2)-°■°^ 
If ~ 10 km s~^ on small sca l es in the dense neutral gas, 
as suggested by lOstriker & Shettvl (1201 ih. this im p lies t s ~ 
0.01, the value propos ed by Krumholz & McKe^ l l2005l) and 
IKrumholz & Taj ( |2007|) and found in Figure |9] In the self- 




Figure 9. Same as Figure [7] but instead of plotting E/tdyn on tlie abscissa we instead plot E/tfi, where tff is the free-fall time in the star-forming clouds 
of a galaxy - see main text for details. We do not include a SgpR — 14^co/ifI plot as the calculation of % requires a gas mass. Symbols are the same as in 
Figures [7] and[8] The left panel shows the results using a bimodal Xqq, while the right panel shows the results using our continuous Xqq. The solid black 
lines show the relation Ssfr = eff Ejnol/iff with eg = 0.01. 



regulation theory of lOstriker & Shettvl ^201 ih . v,. /eg ~ (l/3)/p x 
p,/m,, so that and eg vary together for a given momentum 
feedback level. 

4.3 Observational Constraints on the Model and Xco Values 
for Observed Galaxies 

In order to assess the validity of our parameterisation of Xco 
(Equation 1101). it is worth comparing our models to the existing 
observational constraints in the literature. 

As discussed in § [T] galaxy mergers at low redshift ap- 
pear to have a range of Xco values, from roughly an or- 
der of magnitude below the Galactic mean to comparable to 
the Milky- Way average (Solomon et al.llT997l : lDownes & SolomonI 
[l998; Bryant & Scoville 1999), though on average the Jf -factor 
from local ULIRGs is observed to be below the Galactic mean 
jTacconi et alj|2008h . At higher redshifts, the constraints on Xco 
from inferred mergers (typically submillime tre-selected galaxies ) 
come from either dynamical mass modeling jXacconi et alj|2008h . 



or dust to gas ratio arguments ('Magdis et al.l201ll) . The infeiTed X- 
factors from high-z: SMGs also appear to be lower than the Galactic 
mean by a factor of ~ 5. 

There are relatively fewer constraints on Xco from high-z 
discs. iDaddi et alj l l2010bh estimated the dynamical masses for re- 
solved CO observations of high-2 BzK disc galaxies. After sub- 
tracting off the measured stellar and assumed dark matter masses, 
they were able to derive an Xco factor by relating the remaining 
(presumably H2) mass to the observed CO luminosity. This method 
recovered a mean X-factor ~ 2 x 10 cm" "2/K-kms"\ This IS 
consistent with the calculation of Xco via dust to gas ratio argu- 
ments for a different BzK galaxy bv lMagdis et al. I( l201lh . 

In order to investigate how our inferred X-factors for observed 
galaxies (utilising our model fit) compare to these determinations, 
on the left side of Figure 1 101 we apply Equation [TO] to the observed 
data from Figure [T] and plot the derived Xco for observed galax- 
ies, binning separately for local ULIRGs, inferred z ~ 2 discs, and 
inferred 2 ~ 2 mergers. First, as a consistency check, we exam- 
ine the inferred X-factors employing Equation 1 1 1 from local disc 
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Figure 10. Left: We utilise our model fit for Xnn fEquation flol to infertile Xqq values for observed galaxies compiled bv lGenzel et alJilOlOt) and lPaddi et alj 
i2010al) . Right: CO line intensity, Wco versus Xqq for the observed galaxies from Figure|2](the colour coding of the points is the same as in Figure|7). The 
Xqq values are both the oiiginal literature values (denoted by 'Xqq (bimodal)", and our derived values (denoted by ''Xqq (continuous)"). 



observations, and denote this by the green dash-dot line in Fig- 
ure [To] As expected, the derived X-factors from local discs form 
a relatively tight distribution around the Galactic mean value of 
2 - 4 X lO^^cm^^/K-km s"\ 

From Figure [To] it is evident that observed SMGs have ex- 
tremely low X-factors, with the bulk of them a factor of a few lower 
than the Galactic mean value. Observed ULIRGs show a large pop- 
ulation of galaxies with lower Xqq values, with a few approach- 
ing the Galactic mean. This is reasonably co n sistent with the range 
of values r eported by Solomon et al. l( ll997h ;lD ownes & SolomoiJ 
( 11998 ') and'Bryant & Scovilld ( Il999l) . 

As mentioned, there are far fewer constraints on Xqq from 
high-2 discs, with the only constraints placing the X-factors near 
the Galactic mean value. The values for some of the observed 
galaxies in Figure[TO]are consistent with these determinations. This 
said, there is a peak in our inferred X-factors for high-z discs at 
values in-between present-epoch ULIRG X-factors and the Galac- 
tic mean value. Our models therefore predict that attempts to derive 
Xco for a larger sample of 2 ~ 2 disc galaxies will indeed iden- 
tify some that have X-factors more comparable to local ULIRGs. 
This means that the true H2 gas masses from observations of high- 
z discs and BzK galaxies may be lower than the values reported in 
the literature, though by factors of no more than ~ 2, depending on 
the assumed X-factor 

Finally, on the right side, we plot the CO line intensity, W^co 
versus both the original assumed Xco values for these observed 
galaxies, as well as our derived values. The purpose of this plot 
is to show how different galaxy populations are moving between 
the centre and right panels of Figure |7] Again, by and large, the 
mergers at present epoch and 2 ~ 2 tend to keep derived X-factors 
lower than the Galactic mean, though there are some that fall to 
intermediate values. 



5 DISCUSSION OF OUR MODEL IN THE CONTEXT OF 
OTHER MODELS AND EMPIRICAL RESULTS 

A number of groups have begun to address problems related to 
the CO-H2 conversion factor in a variety of astrophysical environ- 
ments. As such, in order to build a more complete theoretical pic- 
ture of Xco, it is worthwhile understanding how our models fit in 



with respect to these models. The methodologies vary from fully 
numerical to empirical. 

iGlover et al.l ( 1201 Ol) and iGlover & Mac Low! ( 1201 ll) devel- 
oped magnetohydrodynamic models of GMC evolution with time- 
dependent chemistry to follow the formation and destruction of 
both H2 and CO. These models showed that H2 can survive even 
in the presence of a photodissociating radiation field at relatively 
low column densities owing to self-shielding, whereas CO can be 
destroyed more easily. 

These mode ls were expanded by IShettv et al.l ( 1201 Ibl) 
and IShettv et al.l j ioi la) who coupled the simulations of 



IGlover & Mac Lowl 1 201 1.) to large velocity gradient radiative 
transfer calculations. Shetty et al. found that model GMCs with 
mean densities, temperatures and metallicities comparable to 
Galactic clouds result in X-factors comparable to the Galactic 
mean value of ~ 2 — 4 x 10^" cm~'^/K-km s~^. Similarly, these au- 
thours found that low metallicity GMCs result in X-factors larger 
than the Galactic mean. These results are similar to what we have 
found in this study, and have a similar physical reasoning: in low 
metallicity regions in their model GMCs, background UV photons 
photodissociate CO molecules while the H2 self-shields to survive. 

In the large S h2 regime, our Xco results and those of 
IShettvetal.l(l20Ilbll3) differ. Along lines of sight with high-surface 



densitv. lShettv et al. 



( l20Ilah show that, due to line saturation, Xco 



begins to rise (see their Figure 2). While there is quite a bit of dis- 
persion between model galaxies, we generically see a decrease in 
mean Xco at high gas surface densities. The differe nces in Xco 
are a c onsequence of differing underlying models. IShettv et al.l 
( l2011ah model individual GMCs that are resolved on sub-parsec 
scales, whereas our models are of entire galaxies of entire GMCs 
resolved at best to ~ 70 pc. In our simulations, galaxies with higher 
average mean GMC surface densities also typically have larger gas 
temperatures and velocity dispersions as was discussed in § 13.31 
The velocity dispersions in the molecular gas often owe to the tur- 
bulence driven in a galaxy merger, and are larger than would be 
expected if the line widths were solely due to the virial velocity 
of the cloud. Similarly, increased SFRs in the high surface density 
gas drive warmer dust, and consequently gas, temperatures. The in- 
creased temperatures and velocity dispersion allow for larger CO 
line intensities per unit H2 gas mass in this scenario. In contrast. 
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surface density variations in tiie GMCs of IShettv et all ( 1201 lal) do 
not necessarily correlate with larger velocity dispersions and tem- 
peratures in the same way as in a simulated galaxy merger. Thus the 
divergent results in the high Eh2 regime owe to different physical 

processes being mo deled. 

More recently, iFeldmann et al.l ( 1201 1[) have studied the CO- 
H2 conve rsion factor on galax y-w ide scales by coupling the GMC 
models of Glover et al. (2010) and Glove r & Mac Lowl i201lh with 
post-processed cosmological hydrodynamic simulations of galaxy 
evolution. While Feldmann et al. do not have direct information 
regarding the velocity dispersion and line widths (and thus re- 
frain from modeling environments such a s ULIRGs), because they 
utilise the chemical reaction networks of iGlover et al.l ( I2OI0I) and 
iGlover & Mac Lowl ( I2OII1I , they are able to study the effects of 
metallicity on the X-factor with explicit numerical simulations 
(rather than the semi-analytic model fit employed here). Despite 
their somewhat more sophisticated treatment of CO formation, the 
results they obtain are essentially identical to ours: in low metal- 
licity disc galaxies, the X-factor scales as Z' to a power between 
-0.4 and -0.7, consistent with our best fit parameterisation in Equa- 
tionfTOl 

Beyond this, Feldmann et al. relate Xco to Eh2, and find 
variations with scale. On large scales, they find little relationship 
between Xco and Eh2, whereas on smaller scales, Xco increases 
with Eh2. The reasons for this apparent discrepancy between the 
models of Feldmann et al. and us are the s ame in the discussion 
regarding th e models of [ Shettv et alj (1201 laf). 

Finally, lObreschkow & RawUngj ( |2009|) applied a Bayesian 
analysis to literature observational data, aiming to relate the X- 
factor in galaxies to various observables. These authours found a 
similar relationship between Xco and Wco as in this paper. Their 
fit power-law index is ~ —0.31. 



6 SUMMARY 

We have examined the effects of galactic environment on the CO- 
H2 conversion factor in galaxies by coupling simulations of both 
quiescent disc and (merger-induced) starburst galaxy evolution at 
low and high redshift with dust and molecular line radiative transfer 
calculations. While quiescent disc galaxies at low-z with metallici- 
ties around solar tend to have mean Xco values comparable to the 
Galactic mean, we find notable regimes in which Xco may differ 
significantly from this value. In particular: 

(i) In low metallicity galaxies, photodissociation destroys CO 
more easily (due to a lack of dust), whereas the H2 can self-shield 
more effectively. Hence, the amount of H2 represented by CO emis- 
sion rises, and the mean Xco is greater than the Milky Way aver- 
age. 

(ii) In regions of high surface density, the gas temperature and 
velocity dispersion tend to be rather large. The former is due to 
heating of dust by young stars, and energy exchange with gas. The 
latter to the typical origins of high-surface density galaxies: ei- 
ther major mergers, or gravitationally unstable clumps in high-z 
gas-rich discs. In this regime, the velocity-integrated line intensity, 
Woo, rises with respect to Eh2, and there is a net decrease in the 
mean Xco from the Galactic mean value. 

(iii) At high-redshift, gas-rich discs may have gravitationally 
unstable clumps that have moderate velocity dispersions, and high 
gas temperatures (owing to elevated star formation rates compared 
to the Milky Way). These galaxies have X-factors ranging from 
the Galactic mean to a factor of a few lower. Our results for rapidly 



star-forming discs at 2: ~ 2 show that for a given set of physical 
conditions, gas-rich discs at high-2 have comparable X-factors to 
mergers at low-z. Discs and mergers are not inherently different 
with respect to their X-factors. Rather, current local conditions de- 
termine Xco for a particular galaxy. 

These results allow us to develop a fitting formula for Xco in 
terms of gas metallicity and CO line intensity that varies smoothly 
(EquationllOl). Applying this formula to CO detections of galaxies 
at both low and high-z, we find that all versions of the Kennicutt- 
Schmidt star formation relation (Esfr — Emoi, Esfr — Emoi/idyn 
and Esfr — Emoi/iff) niay be described as unimodal with signif- 
icantly reduced scatter compared to literature results that use bi- 
modal Xco values for discs and mergers. 

The results from this work will enable relatively straightfor- 
ward application of an X-factor which varies with galactic envi- 
ronment based simply on two observable parameters (metallicity, 
and velocity-integrated CO line intensity). 
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APPENDIX A: DETAILS OF NUMERICAL MODELING 

In this section, we detail the methodology governing our hydrody- 
namic modeling, calculations regarding the chemical and physical 
state of the ISM, and radiative transfer. 



Al SPH Galaxy Evolution Simulations 

We run GADGET-3 smoothed particle hydrodynamic (SPH) simu- 
lations for both model disc galaxies and galaxy mergers. The aim 
of these calculations is to return the spatial distribution and ages of 
stars, as well as the metal content and gas content as a function of 
location in the model galaxies. 

GADGET-3 is a modified version of the publicly available 
GADGET-2 which employs algorithms for feedback from active 
galactic nuclei (AGN), and better load balancing on parallel pro- 
cessors. For a full description of the underlying algorithms of 
GADGET-2, see lSpringeil ( l2005h . 

The ISM is modeled a s multi-phase, with cold clouds em- 
bedded in a hotter phase jMcKee&Ostrikerl Il977h which is 
practically implemented in the code via hybrid SPH particles 
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dSpringel & Hemauislll2003h . The phases exchange mass via ra- 
diative cooling of the hotter phase, and supernova heating within 
the cold clouds. Within these cold clouds, stars form according 
to a volumetric , Schmi dt ( 1959) power-law relation with index 1.5 
( lKennicuttll998iJ) . The normalisation of this relation is set to match 
the su rface-density Ssfr — Sgas relation as observed in local galax 



ies bv lKennicutll 1 1998a||b|). E xperiments have found that imposing 



a volumetric Schmidt (1959) relation results in a power-law rela- 
tion c omparable to the observed one for the sorts of galaxies studied 
here iMihos & Hemauistlll994al:ISpringelll2"oOol : ICox et ai]|2006bl 



though se elSchave & Dalla Vecchis ( l2008l) for cases when this may 
not apply Ij. 

The gas is initialised as primordial, with all metals forming as 
the simulation progresses. A mass fraction of stars (consistent with 
a Salpeter IMF) are assumed to die instantly in supemovae, and 
enrich the surrounding ISM with metals via an instantaneous re- 
cyclin g approximation using a yield of 0.02 ( Springel & Hernquis3 
I2OO3I ). Hence, gas has a non-zero metallicity from the first genera- 
tion of star formation. We initialise in this manner in order to probe 
a large dynamic range of metallicities in the simulations for both 
disc galaxies and mergers. The typical final metallicities of each 
simulation examined are ~ solar, and we list these in the Table. 

Similarly, supernovae impact the surrounding ISM via energy 
deposition. This pressurisati on of the ISM is implemented via an 
"effec tive" equation of state ( [Robertson et alj|2004l : ISpringel et al.l 
I2OOA . In the models presented here, we assume a relatively modest 
pressurisation of the ISM (qeos = 0.25) for the model z=Q galax- 
ies, and more ext reme (qeos = 1 ) for th e z=2 model galaxies. 
While the tests of [Narayanan et all ( l2011bl) show that o ur results 
are ins ensitive to the choice of equation of state within the lSpringell 
( I2OO5 ) formalism, we are intentional with our choices. The large 
pressurisation of the ISM for the high-redshift models is chosen 
to prevent runaway fragmentation in the dense, gas-rich environ- 
ment (to be discussed shortly). Because our simulations are non- 
cosmological, and do not include gas accretion from the IGM (nor 
d o they have algorithm s for the inclusion of hot gas in the halo as 
m lMosteretal.l20Ilal lbl). extreme star formation in the early phases 
of a galaxy merger would deplete the gas supply prior to the galaxy 
merging. In this limit, galaxy mergers at high-redshift would not 
undergo a starburst. Consequently, we employ a rather stiff EOS. 
We note that while this reduces the amount of fragmentation in 
the ISM ( ISpringe]||2005l) . large ~kpc-scale clumps do still form 
and are dynamically unstable even in t he model disc galaxies at 
z ^ 2, similar to both observations (e.g. Elme green et al ] |2009bl l3. 



(e.g. Bournaud et al. 


2008: Dekel et al. 2009l;lBoumaud et alj201(]|: 


ICeverino et al.i.2010, 


.Bournaud et al.,,201 1 ). 



Black holes are included in the simulations as sink particles 
which accrete according to a Bondi-Lyttleton-Hoyle parameteri- 
sation jBondi & Hovle 1944) with a fixed maximum rate corre- 
sponding to the Eddington limit. AGN feedback is included as 
thermal energy deposited by the central black hole. Specifically, 



^ In iNaravanan et alj i201 Ibh . we explored the effects of modifying this 
SFR power-law index, and found our modeled X-factors were insensi- 
tive to variat ions in the SFR index s o long as the index is > 1. Be- 
yond this, in [Narayanan et alj j201 Ibh . we argued that a Schmidt index 
of unity is unlikely to describe the starburst environments modeled here 
as this choice of power-law inde x does not allow mergers to undergo a 
starbu rst. Similarly, observational iBigiel et al.''2008') and theoretical work 
iKrum holz et al. 2009b; Ostriker & Shetty 2011) suggest that the SFR in- 
dex may be superlinear in high-surface density environments. 



0.5% of the accreted mass energy is reinjected spherically into 
the surrounding ISM which was chosen to allow merger remnants 
to match the present-day A/BH-.A/buisc relation ( Di Matteo et alj 
[2005; Springel et al. 2005; Hopkins et al. 2007a, 2 008b). As dis- 
cussed in the Appendix of [Narayanan et all j20I Ibh . the inclusion 
of AGN feedback does not substantially impact the modeling of 
Xco from galaxies in our simulations. 

The discs are exponential, and initialised according to the 
IMo et alj ( IT99&) model, and embedded in a live cold dark mat- 



ter halo with a Hernquis^ jl990l) density distribution. The mergers 
simply involve discs constructed in this manner. The initial bary- 
onic gas fractions are set to 40% for the z=0 simulations, and 80% 
for the high-redshift models. Because the galaxies consume their 
gas rapidly, the bulk of the high-redshift snapshots which are anal- 
ysed in this study have gas fractions ranging from fg — 0.2 — 0.6, 
comparable to rece nt CO measurements of 2: ~ 1-2 galaxies by 
iDaddi et al.l j20I0bh and|^cconi et al. (2010). 

The halo concentration and virial radius for a halo of a given 
mass is motivated by cosmological TV-body simulati ons, and scaled 
to ma tch th e expected redsh i ft-evo lution following iBuUock et al.l 
( l200lh and [Robertson et"ZI j2006h . The gravitational softening 
lengths are set at 100 h~^pc for baryons and 200 h^^pc for dark 
matter. We simulate a wide range of galaxy baryonic masses, 
merger mass ratios and merger orbits for both z — and z — 2 
galaxie|3- In the Table, we summarise the physical parameters for 
the galaxy evolution models employed for this study, as well as the 
reasoning for the inclusion of that particular model. 

We note that recently the reliability of SPH for cosmological 
modeling has been called into question by comparisons between 
simulations done with GADGET- 3 and others done with the new 
movin g mesh code AREPO deve l oped by[S Pringel (2010). In par- 
ticularJVogelsberger et~al] i20l lh;ISiiacki et a l. (201 1); Kereset^ 



( I2OI ih and lBauer & Springell ^201 ih" have shown that under certain 



conditions, broadly in circumstances where gas in very different 
phases is in close proximity or relative rapid motion, the results 
with the two codes can be very different. In our galaxy-scale simu- 
lations, however, we do not attempt to resolve the different phases 
of the ISM directly, and instead rely on subresolution models to 
describe the gas on small scales, a strategy dating back to early 
work by Hernquist (1989). For this reason, galaxy merger simula- 
tions done with AREPO yield results for the star formation history 
and mass fraction in newly formed stars that are virtually identi- 
cal to those obtained with GADGET-3. Therefore, we believe that 
the results described here are robust with respect to the numerical 
algorithm used to perform the hydrodynamics. 



A2 Dust Radiative Transfer 

We utilise the publicly available SUNRISE dust radiative trans- 
fer simulation package (^Jonsso i^ l2006l : IJonsson et al.ll2006l I2OI0I ; 
lonsson & Primack 2010) for two purposes: to generate the adap- 
tive mesh from the GADGET-3 simulations on which we run the 
radiative transfer, and to calculate dust temperatures in the ISM of 
the model galaxies. 

In constructing the mesh, the physical conditions are projected 
onto a 5'^ base grid spanning 200 kpc. The cells recursively refine 
into 2'^ sub-cells based on the refinement criteria that the relative 



The high-redshift galaxy simulations are actually initialised at z=3. This 
was chosen so that after the ~Gyr typical for the galaxies to reach final 
coalescence when merging, the epoch would correspond roughly to 2 ~ 2. 
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Table Al. Galaxy evolution simulation parameters. Column 1 refers to the model name. 
Column 2 is the baryonic mass in Mq . Column 3 is the merger mass ratio, with "N/A" 
denoting when a galaxy is an isolated disc. Please note that we bin all mergers in the text 
as 1:1, 1:3 and 1:10, though these are approximate. Some l:3-designated mergers may 
be in reality closer to 1 :4, for example. Columns 4-7 refer to the orbit of a merger (with 
"N/A" again referring to isolated discs). Column 8 is the initial baryonic gas fraction of 
the galaxy. Column 9 refers to the redshift of the simulation, and Column 10 is the final 
mass-weighted mean metallicity of the system in units of solar metallicity. 
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density variation of metals Op-m/ < pm. > should be less than 0.1 
and that the l/-band optical depth across the cell is less than unity. 
The maximum refinement level is 11, so that the smallest cells on 
the grid are of order ^ 70 pc across. 

For the dust radiative transfer, the sources of light are stellar 
clusters and the AGN . The stars emit a STARBUR ST99 spectrum 
jLeitherer et al.lll999l ; Ivazauez & Leitherej|2005l) with ages and 
metallicities known from the SPH simulations. The AGN emits 
a luminosity-dependent S EP which is based o n observations of 
unreddened type I quasars jHopkins et al. l l2007hV The normalisa- 
tion of the input spectrum is set by the total luminosity of the black 
hole(s). 

The radiation from stars and the AGN traverses the dusty ISM, 
and is scattered, absorbed, and re-emitted. The evolving dust mass 
is set by assuming a constant dust to metals ratio comparable t o 
that of the Galaxy ijp wek 199 81: 1 Vladildl 1 9981 ; ICalura et al.ll2008h . 
We use the IWeingartner & Draine j200lh" dust grain model (with 



R = 3.15) as updated bv lDraine & H ( |2007|) . The dust and radi- 
ation field are assumed to be in radiative equilibrium, and the dust 
temperatures are calculated iteratively. 

The low-redshift and high-redshift models differ with re- 
gards to the specification of their ISM properties in the SUNRISE 
dust ra diative transfer calculations. Following 'Narayan an et al.l 
( l201Ibh , for the low-redshift models, we assume that young (< 10 
Myr) stellar clusters are embedded in photodissociation regions 
(PDRs) and HII regions for some fraction of their lives. In this 
case, the STARBURST99 spectrum from the stellar clusters is re- 
placed by SEDs der ived from MAPPINGSIII photoionisation models 
( iGroves et alj2008l) . The time-averaged covering fraction of PDRs 
is a free-parameter, and is set to be ro ughly 2-3 Myr . This value 
is motivated in part by simulations of Ijonsson et al.l ^2010 1 who 
showed that this parameter choice when applied to models similar 
to ours produces synthetic SEDs comparable to those observed in 
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the Sp itzer Infrared Nearby Galaxy Survey sample jKennicutt et alj 
I2OO3I) . 

For very massive, gas-rich galaxies, during a merger the stellar 
densities can become so high that stellar clusters begin to overlap. 
At that point it makes more sense to describe the birthclouds that 
the clusters live in as surrounding the entire overlapping (super- 
)stellar cluster structures that form (sometimes as massive as many 
X 10* Mq). The MAPPINGSIII photoionisation models saturate for 
cluster masses of these sizes, and thus become inappropriate to use. 
Therefore, for the high-redshift calculations, we abandon the birth- 
cloud model previously described, and assume the cold ISM is a 
uniform medium with a volume filling factor of unity. While in real- 
ity the ISM may be patchy on these scales, without any information 
as to what this patchiness may be like on resolution scales smaller 
than the SPH smoothing length, we are forced to make the simplest 
assumption possible. That said, this may be a reasonable choice. 
Some evidence exists that in local mergers, a molecular ISM with 
a volume fillin g factor may blanket young stars for the majority o f 
their lives (e.g. lDownes & Solomonll 19981 ; ISakamoto et al.lll999h . 
In any case, t he choice is made prima rily due to simulation code 
constraints. In'Naravanan et all j20I Ibh , we ran test cases examin- 
ing the effects of both varying birthcloud clearing time scales, as 
well as the uniform volume filling fraction model, and found mini- 
mal effects on the derived X-factors. 



A3 The Thermal and Chemical State of the Molecular ISM 

We assume the neutral mass in each cell in the adaptive mesh 
is locked into a spherical isothermal cloud at constant den- 
sity. The surface density of the cloud is that which the simula- 
tion returns, though we impose a floor surface density of 100 
Mq pc~^. This threshold surface density, comparable to the sur- 
face density of most molecular cl ouds in the Local Group (e^. 
Solomon et al. 1987; Rosolowskv e t^l2003l ; Rosolowsky 2005, 
2007 1 ; iBlitz & Rosolowskv. .2006. ; .Blitz et alj [2007. ; .Bolatto et al., 
20081) . is imposed to prevent clouds from having artificially low 



surface densities in large cells in the adaptive mesh. When clouds 
exceed this floor value, we consider them resolved. 

We determine the H2 -HI balance in t hese c louds utilis- 
ing the analytic model o f iKrumholz et al. I ( l2008l. l2009ah and 
iMcKee & Krumholj l2oTol) which treats the balance between pho- 
todissociation of H2 molecules by Lyman- Werner band photons 
against the formation of molecules on grains. The equilibrium 
molecular fraction is given by: 



/h2 



4 1 + 0.25s 



(Al) 



for s < 2 and /h2 = for s > 2. s = ln(l + 0.6x + 
0.01x^)/(0.6rc), where x = 0.76(1 + S.IZ'^'^'^^), and Tc = 
0.066Ecioud/(MQ pc~^) X Z' . Z' is the metallicity divided by the 
solar metallicity. A comparison of this meth od against fully-time 
depen dent chemical reaction networks by Rrumholz & GnedinI 
( I2OI ih finds good correspondence between the methods above 
metallicities of Z' ~ 0.01. While we explore the effects of metal- 
licity on the X-factor in galaxies in these models, we never con- 
sider cases with metallicities lower than this value. 

With the mass and surface density of the cloud known, the 
volumetric density is as well. We scale this density by a factor e'^p^'^ 
to account for the turbulent compression of gas, where numerical 
simulations suggest that 



is a good approximation, where Mip is the 1 dimensional Mach 
number of the turbulence 1 Ostriker et al .l2q0ll;|Padoan & Nordlundl 



l2002l ; lLemaster &Stonell2008l;|Priceet alJl201ll) . The Mach 
ber is calculated by assuming the cloud temperature is 10 K (to 
avoid having to iterate as the temperature we will calculate is de- 
pendent on the density of the gas). 

The velocity dispersion of the gas is determined as the mean 
square sum of the subgrid turbulent velocity dispersion and the re- 
solved nonthermal velocity dispersion. The subgrid dispersion is 
calculated from the external pressure, P — pcciifT"^, and has an im- 
posed ceiling of 10 km s~^, as suggested from galactic si mulations 
which model turbulent energy driving and dissip ation jPib et all 
I2OO6I ; IJoung eTai] |2009| ; lOstriker & Sh"ettvll20Ilh . The resolved 
nonthermal component is calculated from the standard deviation 
in the velocity dispersion in the nearest neighbour cells in the x, y 
and z directions. In cases where the cloud is unresolved, we sim- 
ply assume the GMC is in virial balance, with virial parameter 
Qvir = 5a^i^R/{GM) of order unity so that 



2.2 km s" 



M 



105 M0 



1/4 



(A3) 



where M is the mass of the cloud. 

Finally, we calcul ate the temperature of t he GMC based on 
the model developed bv lKrumholz et al.l tOl Ibl) . In this, we calcu- 
late the heating and cooling processes on the gas, the heating and 
cooling of the dust, and the energy exchange between the two. The 
heating processes of the gas are the grain photoelectric effect, and 
cosmic rays, and the cooling occurs via either CO or CII line emis- 
sion. The dust is heated by the background radiation field, and cools 
thermally. Formally, if we denote heating processes by F, cooling 
by A, and energy exchange via VP, then we have: 



Fpe + FCR - Aline + *gd = 
Fdust — Adust — '^gd = 



(A4) 
(A5) 



For brevity, we refer th e reader to iKrumholz etall j2011bl) and 
[Narayanan et aljfcol l b') for the equations regarding the photoelec- 
tric effect and cosmic ray heating terms. 

The gas cooling is assumed to happen either via CII or CO 
line emission. Following ' Wolfireetal1j20ld ). we approximate the 
fraction of hydrogen where CO is the dominant form of carbon by: 



-4(0.53-0.045111 — 

/CO = ,/h2 X e "h/<: 



-0.0'J7lnZ')/A^ 



(A6) 



which is consistent with the numerical work of lGlover & Mac Low! 

l|). When fco is above 50%, we assume the cooling happens 
via CO line cooling. Otherwise, CII is the dominant coolant. The 
extinction is converted via: 



Av = 



1.87 X 102icm-2 



:Z' 



(A7) 



ln(l + 3Mfo/4) 



(A2) 



(e.g. IWatsonl l20Tlh. We assume the dust is heated primarily by 
background infrared radiation as heating by UV radiation is likely 
to be highly suppressed by extinction. The IR radiation field is 
known from the SUNRISE dust radiative transfer modeling. The es- 
cape of the photons from the GMCs is calculated via th e public 
escape probability code of IKrumholz & ThompsonI ( l2007h . which 
we describe in the subsequent section along with the other molec- 
ular line radiative transfer equations. 
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A4 Molecular Line Radiative Transfer 

Once we know the physical state of the GMCs in our model galax- 
ies, we perform molecular line radiative transfer to calculate the 
CO intensity from the clouds. We first utilise an escape probability 
formalism to determine the escape fraction of photons from indi- 
vidual clouds, and then 3D non-LTE Monte Carlo molecular line 
radiative transfer calculations in order to calculate the transport of 
photons through the galaxies. 

The escape probabiliti es are c alculated using the publicly 
available code described in iKrumh olz & Thompson ( 2007). The 
levels are assumed to be in statistical equilibrium, and are calcu- 
lated by balancing the excitation and deexcitation of CO by colli- 
sions with H2 and He, absorption, stimulated emission and sponta- 
neous emission via the rate equations: 



fu 



(A8) 
(A9) 



where C are the coUisional rates, / the fractional level populations, 
A are the Einstein rate coefficients, and Pui is the escape probabil- 
ity for transition u I. The rate equations are re-arranged as an 
eigenvalue problem and solved with public packages available in 

the GNU Scientific Library. 

P is well-fit by the approximate relation jPrainelbOllh 



where the optical depth is 



gi 16(27r)3/2cr 



1 + 0.5r„ 



QiVH2/i 1 



(AlO) 



(All) 



where Q is the abundance of CO with respect to H2, gi and gu are 
the statistical weights of the levels, 7Vh2 is the column density of 
H2 through the cloud, A„; is the wavelength of the transition, and 
a is the velocity dispersion in the cloud. We iterate Equations |A8]- 
I Al II with standard Newton-Raphson methods until /3 and the pop- 
ulations are known for all levels simultaneously. 

We then propagate the photons which es cape the mode l GMCs 
through the galaxy following the methods of lSernej jl979tF^ . The 
intrinsic line profile function is assumed to be Gaussian in nature 
with width given by the Doppler-width. The model photons are 
emitted isotropically, and can be absorbed by GMCs they encounter 
along the path with a probability determined both by the volume 
filling factor of the GMC in the cell, as well as the absorption line 
profile in the cell (thus accounting for decreasing effective optical 
depths owing to the absorption profile and emission profile shifting 
out of resonance when the line of sight velocity difference between 
the emitting and absorbing cell is large). When all the model pho- 
tons have been emitted, the rate equations are again evaluated, level 
populations updated, and new model photons are emitted. This pro- 
cess is iterated upon until convergence (a fractional difference in 
the level populations of 1 x 10""^) is reached. In practice, this hap- 
pens relatively quickly as velocity gradients in the galaxy render 
the ISM globally optically thin to most photons. 



In lNaravanan et al] ilOl Ibl) . we fully describe the equations employed 
for the galaxy-wide radiative transfer. The interested reader should refer 
to this paper for more detail, and here we summarise the relevant physical 
processes for the sake of brevity. 



APPENDIX B: SIMULATED OBSERVATIONAL 
PROPERTIES OF MODEL GALAXIES 

Because our methodology involves dust radiative transfer, we are 
able to make some comparisons to observed properties of galaxies. 
Here, we present a few comparisons to observations primarily to 
demonstrate that our simulated galaxies serve as reasonable analogs 
to the sorts of galaxies typically observed in large CO surveys. 

At high-redshift, the bulk of the galaxies observed in CO are 
identified according to two selection techniques: the BzK popula- 
tion, and submillimetre-selected galaxies. BzK galaxies are defined 
by their optical colour ratios, and galaxies with blue B ~ z colours, 
and red z — K colours are typically denoted as star-forming BzK 
galaxies due to the presence of a Balmer break. In Figure IbT] we 
examine the optical colours of our model high-z discs. We redshift 
the discs to z=2, and plot the observed {B — z) and (2 — K) colours 
as a compar i son to z ~ 2 BzK galaxies. As originally discussed by 
iDaddi et al.l ( |2004|) . galaxies which fall above the line denoted in 
Figure | B1 1 tend to be disc-li ke in morphology i lDaddi et al.|[2005l : 
iFQrster Schreiber et al.ll200^ . As is evident, our three model high- 
z discs all fall in the same space as observed BzK galaxies. 



[Narayanan et al.l ( l2010bl) and lHavward et al] l l201lh presented 
a merger-driven model for the formation of high-z submillimetre 
galaxies utilising merger-models similar to our high-redshift merg- 
ers. These galaxies were seen to match the mean SEDs of ob- 
served SMGs, as well as the typical CO excitation patterns, line 
widths, and morphologies (Naravanan et a l. 2009). Furthermore, 
when convolving submillimetre-luminous duty cycles with theo- 
retical galaxy merger rates and observed galaxy mass functions, 
iHavward et all (2010) and C. Hayward et al. (in prep.) find that 
this merger-driven model provides a reasonable match to the ob- 
served number counts and redshift distribution of SMGs. In this 
sense, our model high- 2; mergers are likely reasonable analogs for 
the sorts of high-z starbursts observed in C O (e.g. ^Greve et alj 
2005; Tacconi et al. 2006; Bouche et al. 2 0071 ; iTacconi et al.,20oi ; 
feothwell et all201(]i ; lGenzel et al.l2010l) . 

Our simulated mergers at low-redshift are similarly reasonable 
models for present-epoch ULIRGs. In Figure lB2l we show the evo- 
lution of the star formation rate, bolometric luminosity, and IRAS 
25iJ.m/60|i.m flux density r atios for merg er model z0d4e, the fidu- 
cial merger from Naravan an et alj ^201 Ibl ). As is evident, when the 
galaxies merge (around T « 0.85 Gyr), the SFR undergoes a burst 
of comparable magnitude to the most heavily star forming galaxies 
locally. Similarly, the bolometric luminosity rises to the point that 
the galaxy would be selectable as a ULIRG. The IRAS infrared 
colours trans ition from "cool" to "warm" as the starburst and AGN 
heat the gas ^Younger et al.ll2009h . 

Beyond these few examples, a number of other compar- 
isons to observations exist in the literature for this exact same 
set of simulations. For the high-redshift sample, the galaxies have 
been shown to reproduce observed pro perties of Spitzer-selected 
24iJ.m s ources |Naravanan et al.ll20I0al) and properties of bright 
quasars faopkins etalj|2006j|bh . The low-redshift galaxy merg- 
ers have been shown to reproduce kinematic properties of early- 
type galaxies jCox et a l. 2006b), observed X-ray properties o f and 
metallicities in mergers iCox et alj|2006al ; iTorrev et alj 201 ih. th e 
broad-band colours of post- starburst galaxies iSnvder et alj|201 11) . 
and the structural properties o f merging and elliptical galaxies 
( iHoDkins et alj2008clldllal[2009l) . 
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Figure Bl. Synthetic {B — z) and (z — K) colours of model high-z disc 
galaxies. Each red squai'e represents an individual snapshot from models 
z3isob6, z3isob5 and z3isob4. The model galaxies are redshifted to 2 = 2, 
and the colours are observed-frame. Typically, galaxies above the solid line 
are observed to be star-forming discs (Daddi et al. 2004), and are represen- 
tative of the sorts of "normal" (not undergoing a burst of star formation) 
high-z galaxies observed in CO. Our model discs have optical observed- 
frame colours comparable to high-z discs. 
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Figure B2. Evolution of SFR, Lboi, and IRAS colour ratios for an example 
local merger (model z0d4e). Upon merging, the model galaxy undergoes a 
vigorous starburst, is selectable as a ULIRG, and exhibits IRAS 25/60 ^im 
ratios compai'able to many local ULIRGs. Therefore, our simulated mergers 
are reasonable analogs for observed ULIRGs in the local Universe. 

APPENDIX C: ENHANCED COSMIC RAY lONISATION 
RATES IN STARBURSTS 

In this section, we assess the vaUdity of our assumption of a Galac- 
tic cosmic ray ionisation rate. Cosmic rays may affect Xco in 
galaxies by providing a source of heat for the molecular gas. Re- 
calling the methodology presented in this Appendix, we assume a 



Galactic cosmic ray ionisation rate in all model galaxies. It is con- 
ceivable that this assumption may break down in starbursts. Ob- 
servations of M82 by the VERITAS collaboration suggest a link 
between star formation and cosmic rays ( Acciari et al. 2009). Sim- 
ilarly, detections of 7-rays from the Galaxy, LMC, NGC 253 and 
M82 by the FERMI group find a reasonable correlation between 
the 7-ray luminosity, and th e product of the supernovae rate and 
gas mass ( Abdo et aill2010a l. Since 7-rays are the product of cos- 
mic rays with hydrogen atoms, a reasonable scaling may be that 
the cosmic ray ionisation rate is linearly related to the star forma- 
tion rate. Indeed, observations of local st arbursts have suggested 
a pote ntial for cosmic-ray heated gas (e.g. lHailev-Dunsheath et aU 
l2()08l) . 

In order to test this, we have run a series of simulations 
in which we increase the cosmic ray ionisation rate linearly 
with SFR. We a ssume a SFR of 2 MQyr"^ for the Milky Way 
( Ro bitaille & W hitney 2010). We span the range of physical condi- 
tions in our models, running a high-z 1 : 1 merger, high-2: disc, low 
z \: \ merger, and low-z disc simulation. 

In these examples, an increased cosmic ray ionisation rate 
which scales linearly with the SFR does little to change Xco- 
We look at a limiting case to illustrate why this is. In merger- 
induced ULIRGs, where the SFR/cosmic ray flux can be a factor 
of as much as 200 higher than the Galactic value (e.g. Submil- 
limetre galaxies), the minimum gas temperature imposed by cosmic 
ray-heating is ~ 30 — 50 K, depending on the density of the gas 
tPapadop oulos et al.ll201 ih . In these environments, the dust tem- 
perature is typically not very different from this. For example, in 
[Narayanan et alj ( I2OI Ibl) . we found the mass-weighted mean dust 
temperature of a local merger was ~ 70 K. In short, the scenarios in 
which the cosmic ray ionisation rate may be increased have a large 
gas temperature anyway due to high dust temperatures and efficient 
energy exchange between the two at high densities. We therefore 
tentatively conclude that enhanced cosmic ray ionisation rates do 
not affect Xco in starburst environments, though we caution that 
if the cosmic ray heatin g rate rises much fast er than linearly with 
SFR (see, for example, [Bradford et al.ll2003h. then this inference 
may b reak down. Based on the calculations of lPapadopoulos et al.l 
( I2OI ij), we estimate that the cosmic ray heating would affect our re- 
sults significantly only if the ratios of the cosmic ray heating rates 
in mergers and isolated galaxies were at least 10 times greater than 
their ratios of star formation rate. 
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